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ABSTRACT 

We present a spectral analysis of a sample of late L dwarfs. We use our latest model atmo- 
spheres and synthetic spectra and optical and K band spectra to determine effective temperatures. 
We derive effective temperatures of 1400 K to 1700 K for L8 to L6 dwarfs. The analysis demon- 
strates that our recent models that rain out the formed dust completely are applicable to optical 
spectra of late L dwarfs and that more consistent models are needed for intermediate L dwarfs 
and for infrared spectra. We compare the results for the effective temperatures with the temper- 
atures of the onset of Methane formation. Our models predict Methane absorption at 3.3/mi to 
occur at about 400 K higher temperatures than Methane absorption at 2.2/im. This is consistent 
with our data and previous observations which show Methane absorption at 3.3/im but not at 
2.2/im in late L dwarfs. 

Subject headings: stars: low-mass, brown dwarfs - stars: atmospheres - stars: fundamental parameters 
- molecular processes 1 



1. Introduction 

L dwarfs (Martin ct al. 1997, 1998, 1999; Zap- 
atero Osorio et al. 1999; Kirkpatrick et al. 1999, 
2000) are dwarfs cooler than M dwarfs and hotter 
than the so called Methane or T dwarfs. They are 
defined by showing decreasing or no TiO and VO 
absorption in the optical and by not yet showing 
CH4 in the near infrared. There are two classifica- 
tion schemes for L dwarfs, the one by Martin et al. 
(1999) and the one by Kirkpatrick et al. (1999) 
which differ especially for late L dwarfs — which 
are the topic of this work. For the purpose of this 
paper we adopt the scheme by Kirkpatrick et al. 

(1999) . The transition from M dwarfs to L dwarfs 
is observationally very smooth and theoretically 
well explained by gradually decreasing the effec- 
tive temperature. The transition from L dwarfs 
to T dwarfs, however, is not as smooth. The orig- 
inal definition of L dwarfs deliberately left room 
between L dwarfs and cooler objects (Kirkpatrick 
et al. 1999). First estimates on the effective tem- 
peratures of late L dwarfs and on the first known 
Methane dwarf G1229B did not clarify if there is 
a jump in effective temperature and if there is a 
missing population of objects between L dwarfs 
and Methane dwarfs. Recent analyses of L dwarfs 
and to some extend T dwarfs can be found in e.g. 
Allard et al. (1996); Marley et al. (1996); Tsuji 
et al. (1996); Tinney et al. (1998); Burrows et al. 

(2000) ; Basri et al. (2000); Pavlenko et al. (2000); 
Leggett et al. (2001). Those analyses range from 
optical to infrared analyses and they range from 
measuring low to high resolution spectra. How- 
ever, the uncertainty about the effective tempera- 
tures and the transition to the T dwarfs remains. 

One aspect of this transition is the onset of 
Methane absorption. L dwarfs do not show 
Methane absorption in the near infrared, how- 
ever G1229B and even Jupiter do. Recently, Noll 
et al. (2000) reported the detection of Methane at 
3.3fim in late L dwarfs which seems to close the 
gap between L dwarfs and Methane dwarfs. And 
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as discussed in Noll et al. (2000) it even makes the 
naming convention "Methane dwarf" ambiguous. 

Regardless, the effective temperatures associ- 
ated with late L dwarfs and the onset of Methane 
absorption is not clear, yet. Noll et al. (2000) es- 
timated relatively high effective temperatures of 
about 1800 K for late L dwarfs, whereas other 
studies suggest effective temperatures as low as 
1300 K for the latest L dwarfs (Kirkpatrick et al. 
2000). Note that the effective temperature of 
G1229B is estimated to be 900-1000 K (Allard 
ct al. 1996; Marley et al. 1996), but that it is 
not the warmest T dwarf since earlier ones (e.g. 
Leggett et al. 2000; Burgasser et al. 2000; Bur- 
gasser et al. 2001) are known to exist. This work 
aims to clarify the effective temperatures of late 
L dwarfs and the onset of Methane absorption. 

In the next two sections we describe the data 
and the models used in this work, respectively. In 
section 4 we present the analysis which consists 
of an optical spectral analysis, an infrared spec- 
tral analysis and an absolute flux analysis for the 
objects with known parallaxes. We discuss and 
conclude this work in section 5. 

2. The data 

The optical spectra for the stars of this sample 
are taken from Kirkpatrick et al. (1999) and Kirk- 
patrick et al. (2000). The K-band spectra used in 
this work were obtained at the Kitt Peak 4-meter 
telescope using the CRSP spectrograph. The spec- 
tra were reduced using standard techniques. The 
atmospheric water absorption was removed using 
nearby A type stars. See also Gizis et al. (2001) for 
details of the observations. The sample constitutes 
of all late (> L6) L dwarfs that have both optical 
and K band spectra. The criterion L6 was physi- 
cally relatively arbitrary. It was a rather technical 
criterion. We wanted to comfortably include the 
spectral range used in Noll et al. (2000) and we 
wanted to have a spectral range for which we have 
valid models. See Allard et al. (2001); Schweitzer 
ct al. (2001) and sections 3 and 5 for a discus- 
sion of the validity of the applied models. Es- 
pecially the respective sections of this paper will 
make this work also a case study for the validity 
of the AMES-Cond models. 

Some of the stars also have known trigonomet- 
ric parallaxes which are collected and tabulated 
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in Kirkpatrick et al. (2000). All objects and their 
properties are listed in Table 1. 

3. The models 

The models used in this work are the so called 
AMES-Cond and AMES-Dusty models described 
in Allard et al. (2001). A detailed description can 
be found therein. The key features relevant for 
this study are the dust treatment and the latest 
line lists which are summarized in the next few 
paragraphs. 

In the AMES-Cond models the dust is treated 
in the limiting case where the dust has formed 
in equilibrium and has rained out completely out 
of the atmosphere. This effectively removes the 
fraction of elements from the atmosphere that are 
bound in dust grains; see Allard et al. (2001) for 
a discussion. The AMES-Dusty models treat the 
dust in the other limiting case where the dust stays 
in the layer in which it forms and does not rain out 
at all. This effectively introduces a dust opacity 
due to the full dust content. To calculate the dust 
content we use the Gibbs free energy of formation 
from Sharp & Huebner (1990) and calculate the 
equilibrium concentration of all dust species avail- 
able. 

In earlier analyses, the models have been ap- 
plied to low resolution optical and infrared spec- 
tra of M and L dwarfs (Leggett et al. 2001) and 
to optical low and high resolution spectra of early 
L dwarfs (Schweitzer et al. 2001). Earlier mod- 
els have been applied to high resolution spectra 
of L dwarfs as well (Basri et al. 2000). The re- 
sults from these works indicate that the assump- 
tion of complete rain out is valid for late L dwarfs. 
Therefore, we will primarily apply the AMES- 
Cond models in our analysis of the late L dwarfs. 
However, to further investigate the condition of 
complete rain out we also used infrared spectra 
of AMES-Dusty models in sections 4.2 and 4.4 as 
comparison. 

The important opacity sources included in this 
work are the TiO(Schwenke 1998), FeH(Phillips 
& Davis 1993), CrH(Frcedman, private comm., 
see also Allard et al. 2001), H 2 0(Partridge & 
Schwenke 1997) and CH 4 (Husson et al. 1992) 
line lists. This means that all important opacity 
sources for (late) L dwarfs are treated with line 
lists and no longer with JOLA (Just Overlapping 



Line Approximation) opacities. However, we note 
that there are still uncertainties in the data for the 
H 2 line list (Allard et al. 2000). 

The spectra used for this analysis are calculated 
at a resolution of 2 A. They are based on fully 
converged model structures that include depth de- 
pended line profiles and include all the opacities 
for all iterations (see e.g. Hauschildt et al. 1999; 
Allard et al. 2001, for details of the calculations). 

4. Analysis 

4.1. Optical spectra 

We used a \ 2 fitting technique to determine 
the best fitting model spectra. It is described in 
Leggett et al. (2001) and has also been used in 
Schweitzer et al. (2001). As input we used a grid of 
AMES-Cond model spectra from T off =1000 K to 
T off =1700 K in steps of 100 K and from log(.g)=4.0 
to log(g)=6.0 in steps of 0.5 to determine which 
model results in the smallest % 2 value. Since the 
observed spectra have a resolution of significantly 
less than R=1000 the standard resolution of 2 A 
of our model spectra is sufficient for this analy- 
sis. However, the resolution of the observations is 
not high enough to analyze individual atomic lines. 
Therefore, the x 2 method used in this work mea- 
sures the goodness of the fit to the overall spec- 
trum and weighting individual lines would not be 
allowed by the quality of the available data. We 
employ error bars for the measured values of 100 K 
in T e ffand 0.5 dex in log(g). This reflects primarily 
the grid spacing. However, the three best fitting 
models do lie within these error bars. These er- 
ror bars do not take into account any systematic 
errors due to uncertainties of the models. 

The best fits to the optical spectra are shown in 
Figs.l to 7 and the parameters for the best fitting 
models are summarized in Tab.l. As can be seen 
the models fit very well to the observations. The 
FeH and CrH bands at 8600 A and the H 2 bands 
at 9300 A are well reproduced as well as the overall 
shape which is determined by the K I doublet at 
7600 A and the Na D doublet. Also the Rb I 
and Cs I are reproduced within the limits of the 
resolution. 

The determined effective temperatures corre- 
late with spectral type, as expected. The re- 
sults might suggest a jump in effective tempera- 
ture from L7 to L8 by 300K. However, this sam- 
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pic is too small to draw such conclusions, espe- 
cially within the error bars of the fitting process. 
The determined gravities are log(g)=5.0 or larger 
which means the objects are well contracted and 
not young anymore. This is expected for a field 
sample like this 2MASS sample and very similar to 
the results for log(g) obtained in Schweitzer et al. 
(2001). 

4.2. Infrared spectra 

We performed the same type of \ 2 fitting tech- 
nique for the infrared spectra as we did for the 
optical spectra. The model grid for the infrared 
analysis had to encompass effective temperatures 
from 1300 K to 1900 K to obtain good fits. The re- 
sulting fits are shown in Figs. 8 to 14 and the best 
fitting parameters are listed in Tab. 1. For the 
K band spectra the resulting models were signifi- 
cantly hotter than the best fitting models for the 
optical spectra. However, results from previous 
(mostly optical) studies showed that AMES-Cond 
models are appropriate only for effective tempera- 
tures below about 1700 K, whereas our fitting pro- 
cess produced significantly hotter effective temper- 
atures. Therefore, we also tried to fit the AMES- 
Dusty models to the K band spectra, but the re- 
sulting temperatures were about the same and the 
quality of the fits did not improve significantly. 
As a visual inspection showed, the shape of the 
K band AMES-Cond spectra is not very temper- 
ature sensitive and mostly dominated by CO and 
H2O absorption for models with T e ff>1800 K. Be- 
low 1800 K, CH 4 starts showing in the synthetic 
spectra which is further investigated in section 4.4. 

4.3. Absolute flux comparisons 

To clarify on the different effective tempera- 
tures we calculated absolute flux values for both 
synthetic spectra and synthetic magnitudes and 
compared them with observed data. The calcula- 
tion of synthetic absolute flux values hinges on the 
knowledge of the radius of the object. However, 
according to Chabrier et al. (2000), the radius of 
a very low mass object is 0.1 ±0.03 R© for an age 
older than 300Myrs and a mass between 0.01 M 
and 0.07 M©. These are very conservative limits 
which are easily fulfilled for a random field sample 
like the presented 2MASS sample. Therefore, we 
used a radius of 0.1 R© (« 1 Rj up iter) to calculate 
absolute flux values. 



Such a large error in radius will introduce an 
even larger error in flux. However, the reason 
for doing this comparison is to compare flux val- 
ues from the "cool" (1400K) models and from the 
"warm" (1700K) models. At these low effective 
temperatures the flux differences from the effective 
temperatures are much larger than the flux differ- 
ences from the radius uncertainty. Therefore, we 
can use the absolute flux values to decide between 
the results for the effective temperatures from the 
optical fits and from the infrared fits. 

Table 1 lists the observed K5 magnitudes and 
K5 magnitudes calculated from the measured 
trigonometric parallax the synthetic spectra with 
the best fitting parameters and the 2MASS Kg 
response function. As can be seen, using the pa- 
rameters from the optical low resolution spectra 
fitting give excellent agreement with observations. 
The synthetic magnitudes using the results from 
the infrared analysis are significantly brighter than 
the observed magnitudes. 

The K s magnitude of 2MASSs J0850359+105716 
is within the error bars of both the 1700 K model 
and the 1900 K model synthetic magnitudes and 
cannot be used to distinguish between the warm 
and cold models. For the two L8 dwarfs with 
parallaxes it would alternatively require an as- 
sumed radius of about 0.15 R© to favour the syn- 
thetic magnitudes from the infrared fitting over 
the synthetic magnitudes from the optical fitting. 
However, according to Chabrier et al. (2000) this 
would correspond to objects younger than 50 Myrs 
and less massive than 0.02 M© which have match- 
ing effective temperatures. This is extremely un- 
likely for a field sample like this one and would 
be in contradiction to recent age estimates for 
2MASSW J1523226+301456 (G1584C) of at least 
lGyr (Kirkpatrick et al. 2001). 

Figure 15 shows 2MASSW J1632291+190441 
and 2MASSW J1523226+301456 and the respec- 
tive best fitting models in astrophysical flux. The 
only scaling applied to the synthetic spectra are 
the distance and an assumed radius of 0.1 R©. 
As can be seen, the flux agreement is excellent. 
2MASSs J0850359+105716 is not shown since the 
best fitting model (AMES-Cond with T eff =1700K, 
log(<7)=5.5) scaled to the distance and 0.1 R© has 
almost a factor two more flux than the obser- 
vation. This is within the expected error range 
due to the radius uncertainty. Also, T e fj=1700K 
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is already close to the parameter region between 
the valid regimes of the AMES-Cond and AMES- 
Dusty models. In particular, AMES-Cond models 
(when leaving their applicable parameter range) 
have too much flux in the optical due to the lack 
of dust opacity. Again, using spectra with the pa- 
rameters obtained from the infrared fitting results 
would be brighter by more than a factor of two 
(aside from the different spectral shape). 

Since the infrared spectra are not absolute flux 
calibrated and since there are no optical magni- 
tudes available we could not do a corresponding- 
analysis for those quantities. 

4.4. Methane behavior 

As mentioned above, the AMES-Cond models 
start showing CH 4 at 2.2^m for T cff <1800 K. This 
is the main reason why the K band fitting only 
yields effective temperatures above 1800 K. The 
AMES-Dusty models do not show CH4 at 2.2/rni; 
not even at T e ff=1500 K. Below that effective tem- 
perature, the AMES-Dusty models are no longer 
appropriate (see e.g Schweitzer et al. 2001; Allard 
et al. 2001) and their spectra are unrealistic. 

However, late L dwarfs (T o ff«1500 K according 
to our optical and flux analysis) do not show CH4 
at 2.2^m. Yet, Noll et al. (2000) have recently 
reported the detection of CH4 at 3.3/xm for late 
L dwarfs. 

Therefore, we also investigated the 3.3/^m spec- 
tral region and found the onset of CH4 to oc- 
cur at an effective temperature of about 2200 K. 
This is about 400 K hotter than the onset of the 
2.2^m band. The absorption onset at 3.3/zm as 
well the absorption onset at 2.2^m are demon- 
strated in Fig. 16. As can be seen, the CH4 
feature is visible at 2.2/zm for T e gf<1800K and at 
3.3^m for T cfr <2200K. The figure shows models 
with log(g)=5.5, however, the result is the same 
for 4.5<log(g)<6.0. Note, that our models are 
fully self consistent. Unlike the study in Noll et al. 
(2000) our spectra are a direct output from the 
model calculations. 

The previous sections suggest that there is a 
systematic offset in the theoretical K band (and 
most likely other infrared band pass) spectra. 
They seem to show a typical late L dwarf shape at 
about 300-400 K too high temperatures. There- 
fore, the CH4 onset temperature of 2200 K at 



3.3/xm has to be taken with care. We are much 
more confident in the differential result than the 
absolute result. I.e. we are confident that the CH4 
absorption at 3.3/im occurs at about 400 K higher 
effective temperatures than the CH4 absorption at 
2.2^tm. This is consistent with the 3.3/xm band be- 
ing a fundamental band and the 2.2/im band being 
an overtone band which only appears when suffi- 
cient amounts of CH4 have formed, i.e. at cooler 
temperatures. 

5. Summary and discussion 

We have performed a multi- wavelength analysis 
of late L dwarfs. From optical low resolution fits 
we determined effective temperatures for L6 to L8 
dwarfs between 1700 K and 1400 K. From K band 
spectral analysis we determined effective temper- 
atures between 1900 K and 1700 K. Comparing 
absolute flux values (K band magnitudes and ab- 
solute flux calibrated optical spectra), we find ex- 
cellent agreement with the cooler effective tem- 
peratures. The mismatch for the K band spectra 
means that the shape looks too hot, but the flux is 
correct. Most prominently, the AMES-Cond mod- 
els predict CH4 absorption at 2.2/xm which is not 
observed. 

This can be caused by several factors or a com- 
bination of them. For one, the dust is treated only 
in a limiting case. A more correct treatment which 
accounts for gravitational settling will change the 
temperature structure and opacity structure of 
the atmosphere and hence parts of the spectrum. 
As already discussed in Allard et al. (2001) and 
Schweitzer et al. (2001), the treatment of dust has 
a significant impact on the temperature structure. 
This together with the improved dust treatment 
will also change how the elements are distributed 
over the compounds and this results in different 
concentrations of CO and CH4 (the most promi- 
nent opacities in the K band) . Secondly, there are 
known uncertainties in the H2O line list (Allard 
et al. 2000) . The current water opacity distributes 
the flux in the J to K bandpasses incorrectly. 

This has interesting consequences on the CH 4 
bands. We found that the onset of the band at 
3.3^m appears at about 400 K higher tempera- 
tures than the onset of the 2.2/xm band. The num- 
bers for the respective effective temperatures have 
been derived from the AMES-Cond models. Pre- 
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vious studies (Schweitzer et al. 2001; Leggett et al. 
2001) have only confirmed their validity for opti- 
cal spectra of late L dwarfs. It appears from this 
study that the K band shape and the onset of CH4 
seems to be reproduced at too high temperatures. 
As already discussed above, this sort of effect is 
consistent with an unrealistic temperature struc- 
ture which necessarily arises from treating dust 
only in a limiting case. 

The AMES-Dusty models seemed to reproduce 
realistic spectra in the optical and near infrared 
when applied to appropriate objects(Leggett et al. 
2001; Schweitzer et al. 2001). However, as noted 
in this study, they do not predict any CH4 absorp- 
tion. The future models we are currently working 
on (Allard et al., in preparation) will treat the 
dust formation and settling in a consistent fash- 
ion. These models will then have more realistic 
temperature structures and reproduce the spec- 
tral features more accurately. On the other hand, 
the excellent fits of the AMES-Cond models in the 
optical already suggest that the outer atmosphere 
is dust free. 

From the range of previous estimates on the ef- 
fective temperatures of late L dwarfs we derived 
relatively low effective temperatures. This leaves 
little room for objects between classes L and T. 
Although the optical spectra look significantly dif- 
ferent, our study found strong theoretical evidence 
that CH4 absorption already takes place at late L 
type supporting the results from Noll et al. (2000). 
More high signal to noise spectra of L dwarfs in 
the L band and K band are needed to underpin 
these recent results. 
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Fig. 1.— Fits (grey line) to 2M0850+105 (dark 
line). See Tab. 1 for parameters. Telluric features 
have not been removed. See Kirkpatrick et al. 
(1999) for exact locations of telluric bands. 
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Fig. 3.— Fits (grey line) to 2M0825+2115 (dark 
line). See Tab. 1 for parameters. Telluric features 
have not been removed. See Kirkpatrick et al. 
(1999) for exact locations of telluric bands. 
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Fig. 2.— Fits (grey line) to 2M0920+351 (dark 
line). See Tab. 1 for parameters. Telluric features 
have not been removed. See Kirkpatrick et al. 
(1999) for exact locations of telluric bands. 




7000 8000 9000 10000 

Wavelength in A 

Fig. 4.— Fits (grey line) to 2M0929+342 (dark 
line). See Tab. 1 for parameters. Telluric features 
have not been removed. See Kirkpatrick et al. 
(1999) for exact locations of telluric bands. 
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Table 1 

Known and derived parameters for the objects in the sample. 



Optical fitting Ks band fitting 



Name 


Spectral Type 


K S 


Wtrig 




log(9) 


Ks,Sy„ 


T c ff 


log(s) 


K S ,Sy„ 


2MASSs J0850359+105716 


L6 


14.46±0.07 


30.6 


1700 


5.5 


14.49 


1900 


6.0 


14.09 


2MASSW J0920122+351742 


L6.5 


13.93±0.08 




1700 


6.0 




1900 


6.0 




2MASSI J0825196+211552 


L7.5 


13.05±0.04 




1700 


6.0 




1800 


6.0 




2MASSW J0929336+342952 


L8 


14.62±0.11 




1400 


5.0 




1800 


6.0 




2MASSI J0328426+230205 


L8 


14.84±0.13 




1400 


6.0 




1700 


6.0 




2MASSW J1632291+190441 


L8 


13.98±0.05 


16.8 


1400 


5.5 


14.06 


1700 


6.0 


13.28 


2MASSW J1523226+301456 b 


L8 


14.24±0.07 


18.8 


1400 


5.5 


14.30 


1800 


6.0 


13.26 



"from Kirkpatrick et al. 2000 
b Gl 584C 

c The error bar on the synthetic magnitudes is about 0.5 magnitudes from the radius uncertainty of 0.03 R© (see text for details) 
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Fig. 5.— Fits (grey line) to 2M0328+230 (dark 
line). See Tab. 1 for parameters. Telluric features 
have not been removed. See Kirkpatrick et al. 
(1999) for exact locations of telluric bands. 
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7000 8000 9000 10000 

Wavelength in A 

Fig. 7.— Fits (grey line) to 2M1523+301 (dark 
line). See Tab. 1 for parameters. Telluric features 
have not been removed. See Kirkpatrick et al. 
(1999) for exact locations of telluric bands. 



Fig. 6.— Fits (grey line) to 2M1632+190 (dark 
line). See Tab. 1 for parameters. Telluric features 
have not been removed. See Kirkpatrick et al. 
(1999) for exact locations of telluric bands. 
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Fig. 8.— Fits (grey line) to 2M0850+105 (dark 
line). See Tab. 1 for parameters. 
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Fig. 9.— Fits (grey line) to 2M0920+351 (dark 
line). See Tab. 1 for parameters. 
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Fig. 10.— Fits (grey line) to 2M0825+2115 (dark 
line). See Tab. 1 for parameters. 
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Fig. 11.— Fits (grey line) to 2M0929+342 (dark 
line). See Tab. 1 for parameters. 
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Fig. 12.— Fits (grey line) to 2M0328+230 (dark 
line). See Tab. 1 for parameters. 
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Fig. 13.— Fits (grey line) to 2M1632+190 (dark 
line). See Tab. 1 for parameters. 
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Fig. 14.— Fits (grey line) to 2M1523+301 (dark 
line). See Tab. 1 for parameters. 
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Fig. 15. — Models (grey line) and observed 
spectra (dark line) for 2M1632+190 (left) and 
2M1523+301 (right) in absolute flux units. The 
only scaling applied to the models are the paral- 
lax and a radius of 0.1 R Q . 




Fig. 16. — The onset of Methane absorption at 
2.2^tm (left) and at 3.3/zm (right). All models are 
for log(<?)=5.5. T ffis from top to bottom 1900K, 
1800K, 1700K, 1600K (left) and 2300K, 2200K, 
2100K, 2000K (right). The models have been con- 
volved to yield a resolution of R=1000. 
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